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The profile of a Bennett hole induced by laser field in ionic distribution in collisional plasma 
is calculated. Influence of Chandrasekhar's dependence of coefficients of velocity space transport 
on the profile is included into the calculation for the first time. It is found that the hole narrows 
down as the field detuning frequency increases. Physical cause of the effect is the falling dependence 
of Coulomb collision frequency on the ionic velocity. Estimations show that the effect is quite 

• observable under conditions of high-current gas-discharge plasma. 
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The nonlinear spectroscopic technique as a plasma diagnostic tool is a subject of considerable recent interest 
because of its novelty and potential in basic and applied plasma physics. Of particular interest are the spectroscopic 
effects of Coulomb scattering |jj , since it differs radically in appearance from binary collisions in gases [^|,|| . In linear 
absorption spectra Doppler width usually exceeds collisional broadening, so it masks effects of the scattering. However, 
the velocity distribution of particles interacting with a monochromatic field near the resonance is changed significantly 
by the scattering. If the frequency ui of the field is near-resonant with Bohr's frequency u) mn = (E m — E n )/h between 
levels m and n of the atoms with velocity v that satisfies the condition kv = Lu — uj mn = Q, then the induced transitions 
m <-> n occur (k denotes the wavevector). As a result, the narrow peak or dip (Bennett hole [Q) arises in the velocity 
■ distribution of level populations. The steady-state profiles of holes are controlled by both the homogeneous width of 
transition m <-> n and the scattering processes with a change in velocity. Measuring the width, shift and asymmetry 

• of nonequilibrium structures in the distribution with the help of probe-field technique one can extract information 
Cj^ 1 about the scattering. 

Previous calculations of Coulomb broadening were based on the assumption of the constancy of the diffusion 
coefficient in the velocity space. This assumption is applicable when the detuning frequency Q is less than the Doppler 
width kvx, where vt — \f 2Ti / m is the thermal velocity of ions and m is the ion mass. The constant diffusion model 
was enough to interpret Lamb dip measurement in argon laser plasma or another with near- resonance pumping. In 
contrast, farther away from the resonance, at Q > kvx, one must take Chandrasekhar's velocity dependence of the 
diffusion tensor that complicates the evaluation. One way around this problem was shown in linear theory of Dickc 
narrowing in the ionic spectrum, which had been built up for low intensity of the incident wave M . 
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In the present Letter we study a first-order nonlinear problem, the influence of collisions in plasma on the profile 



of the Bennett hole. The calculation presented below for nonlinear structure in the velocity distribution is performed 
in the basis of density matrix formalism ||J7]]. First we derive the formula for the distribution, then analyze it under 
some simple but realistic conditions, and at the end estimate whether the effect obtained be observable in experiment. 
In particular, we show that the width of hole decreases with the detuning of the electromagnetic field from the 
resonance. By observing the effect, one could directly measure the diffusion tensor. It would be helpful to confirm 
C$ ' our understanding and offer scope for new diagnostic methods. 

The diffusion coefficient as a function of the velocity was measured for Q-machine plasma at low density by Bowles 
ct al. [|| . The effect of transient spread of a Bennett hole in the distribution of tagged ion population in a metastable 
state was exploited. However, that effect cannot be used for denser plasma, because the diffusion time becomes 
smaller and the measurement requires a faster registration technique. In this paper we shall examine the stationary 
hole shape. Analysis of its width is useful for diagnostics at high as well as at low density. 

To develop the theory of nonequilibrium structures in the velocity distribution, let us describe the two-level subsys- 
tem of a probe ion in an ideal nondegenerate plasma by the spectroscopic density matrix pij (r , v, t) . This matrix over 
internal states i, j is at the same time the Wigner function of position r and velocity v. It has been shown J^,|^] that ele- 
ments of the matrix satisfy the quantum kinetic equation with the classical Landau collision term. This equation makes 
it possible to analyze resonant transitions between levels in the field of the traveling wave E = Eq exp(— iu>t + ikr) + 
c.c. 

Denoting the relaxation constants of level populations as Tj and that of the coherence as T, write steady - state 
equations for off-diagonal p m „(r, v, t) = p(v) exp(—iQt + ikr) and diagonal pjj = pj elements of the density matrix 



r - ift + ikv) p - vVp = -iGAN, T 3 p 3 ■- vV Pj - \ 3 {v) = T^(iG*p), (1) 
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where the upper/lower sign corresponds to upper/lower level. Here AN = p m — p n , G — Eod mn /2h, d mn is the matrix 
element of the dipole moment, 
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V ~ 3m 2 4 (2) 

is the Coulomb collision frequency, Z a e is the charge of a probe ion, L is the Coulomb logarithm. We assume that 
there are only singly charged ions in plasm; iVj is their concentration. The excitation of ionic levels j — m,n usually 
occurs from the ground or metastable state, so we are justified in assuming the shape of the excitation function Xj(v) 
to be Maxwellian 

\ 3 (v) = Q 3 W{vl ^ = ^«ap(-|)> (3) 

where Qj is the excitation rate of level j . 

For plasma in a state close to the thermodynamic equilibrium with ionic temperature Tj, the distribution of buffer 
particles has the Maxwellian shape. Then it is possible to write an explicit expression for collision operator V involving 
the dynamic friction and the diffusion in the velocity space 

where £ Q = v a /vT is the a - component of the dimensionless velocity, the differential operator d/d£ a acts on the 
right. Functions and <I>t r (£), occurring in the expression for the diffusion tensor 

^ = + <M0 (Safi - ^) , (5) 

can be written as single integrals 

1 l 
*,(0 = 3 J \ 2 e-^ 2 dX, = | f(l - X 2 )e~^ 2 dX. (6) 
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The velocity distribution of populations pj of both levels j = m,n needs to be found. To simplify analysis, we 
assume the electromagnetic field to be weak \G\ 2 <C ITj and will neglect the effects of strong saturation. The right 
sides of Eqs. (0) can therefore be regarded as perturbation. The formal solution of in terms of operator exponents 
is 

pf - y^W(v), pW = 0; AN(v) = p° m (v) - (v) = N mn W{v); 

p {1) (v) = -iG I dtexp(-At + vVt)AN(v), A = T-in + ikv, (7) 
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(v) = T / dt exp(-r i i + vVt)2$t (iG*p {1) («)) . (8) 



It remains to write out the explicit solution and analyze its limiting cases. For this purpose we must calculate 
the commutators of operators A and V and get, in the general case, the infinite series. Fortunately, the relation 
[[[V, A], A], A] =0 allows us to break the series within the first order in v. 
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P y>{v) = ^2\G\ 2 AN(v)U J dt J drexp {-T jT -(T-iCl + ikv z )t) x 



V - (-<S> zz (£)k 2 v 2 T (t 2 T + t 3 /3) + 2ikv z ^(0(2tr + < 2 )) + £V)] , (9) 



x exp 

where the axis z is chosen along the wavevector k. Function 
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<M0 = (MO - MO) | + MO (io) 

is the zz - component of the diffusion tensor. This function depends on two variables £ = |£| and £ z . The integrand in 
(B), obtained by the expansion, is inappropriate at long times t, r. Consequently, for application of this expansion it 
is necessary to have the main contribution to integral (|9|) in the domain where the quadratic terms in v remain small. 
Since \d<fr a p/ 9£ 7 | < |$ a( a|, estimation of those terms gives the following applicability conditions of the expansion (||) 

« = *=£<1. (ID 

The latter inequality means the Doppler limit; the former one coincides with the condition of smallncss of the diffusion 
broadening compared to the Doppler width. The collision frequency in high-current gas-discharge plasma is about 
v ~ 10 5 - 10 7 s^ 1 (see 0). It is by a few orders of magnitude less than the Doppler width kvr ~ 10 10 s _1 and, as a 
rule, less than the relaxation rates T ~ 10 8 — 10 9 s -1 , Tj ~ 10 6 — 10 9 s _1 . That is why we may use small parameters 
a, (3 and it is enough to hold only terms linear in v. 

Taking an integral over r we obtain the correction to the distribution function as a single integral over t 



exp 



-At - $ zz (£)i/k 2 v%,t 3 /6 + ivkv z $i(£)t 2 



p\ L >(v) = t2\G\ 2 AN(vW / L - — — — -±dt. (12) 

3 J T j +$ zz ((\vk 2 v 2 t 2 /2-2ivkv z $t(Z)t 

The integrand in ([l2]) regarded as a function of t at fl = kv z is maximal at t = and decreases with t. The shape 
of the distribution depends on a few parameters. Considering conditions ( [ll"l) and the inequality Tj < T we can set 
apart two limiting cases of parameters: 

vk 2 v 2 T < T 2 T 2 (a < (3 2 ), (13) 



T% < vk 2 v 2 T (a > [3 2 ). 



(14) 



For the domain (|13|) integral ( |12| ) is accumulated at i < Ti = 1/r. The diffusion and friction cause the Bennett hole 
to somewhat change its width and position. For the limit ( |l4| ) integral ( |l2|) is gathered at t < T 2 = (a/2)~ 1 / 2 /kvx- 
Diffusion and friction effects control the width and shift of the hole. The width is defined by the maximum among the 
homogeneous width Avh = T/k and diffusion width Avjj = vt{v /T 'j) 1 ^ 2 ■ Thus, the scattering has stronger impact 



on shape of the hole, while the relaxation constants of selected levels satisfy inequality (14). 

Let us find the width as a function of detuning in this limit e(0 = 2(3 2 / 'a<I> Z2 (0 <C 1. It is convenient to write the 
expression for the nonlinear correction to the Maxwellian distribution as a convolution of two functions 



pf\v) = W(v) f f(x-x') gj (x')dx', 



= kv, — f2, 



(15) 



where 



f(x) = 5ft 



dt 
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M(D=e(0+^ 2 -e(0: 



(r + ix)t - § zz {()vk 2 v^t /6 + ivkv z $i(£)t 



(17) 



In the limit e <C 1, the width of function /(#) is much less than that of gj(x). Accordingly, the central part of the 
distribution is given by expression (|l6| ) pj ~ 9j(^) at |x| < T/yfe. Meanwhile function gj(x) drops exponentially 

with so the asymptotic of the distribution p^\v z ) at \x\ T/y/e is given by a narrow but slowly falling function 
f(x) and has the Lorentzian shape. 
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The most interesting feature of the distribution obtained is the narrowing of its central part as the detuning S7 
increases, due to the falling velocity dependence of the collision frequency. One can estimate the half-width by 
substituting the extreme value of the longitudinal velocity v z = 0,/k into all components of the diffusion tensor: 




This quantity is proportional to the square root of the zz- component of the diffusion tensor, ft falls off with detuning 
S7. The effect is absent in the model of the constant collision frequency fl^JTol- 

To describe this effect quantitatively we calculate the ft- dependence of x 2 averaged over the distribution (|lE|). 
Keeping only the terms linear in parameters a and (3 we obtain 

F s (n) = J d\ x 2 pf(ti) ~ T 2|G ^'"" I r + jdt$ zz (t) exp (-(r - ifi)t) 1 , (18) 



1/V2 

$zz(t) = J d 3 vW(v)$ zz (g)exp(-i£ z t) = J 6y 2 dyexp 

o 

The first term on the right-hand side of Jl8| ) is independent of S7 and corresponds to the contribution of the Lorentzian 
wings of the distribution. The second term describes the collisional broadening and decreases with 

Under the approximations the line shape is Voigt's contour. The dependence of the variance D = ((x — (x)) 2 ) on 
detuning SI consists of two terms: the contribution of the wings, which increases with detuning, and the decreasing 
contribution of the diffusion broadening. At small detuning |fi| <C kvx we have the explicit expression 

Hwr ( f SI \ 2 \ y{kv T f fflW,. 

) + ^-[y- 3 -{i^) 1 15 -tJJ- 

Summarizing the results, we conclude that the Bennett hole narrows down as detuning increases. The physical 
reason of the effect is that the width of the hole is not affected by the homogeneous width Avh, but depends on 
the change in velocity due to diffusion Av& = VTy/^Tj, while Avjj 3> Avh- The excited state lifetime Tj = TJ 1 
determines the stationary width of the hole. While the detuning SI of the held remains inside the Doppler contour, 
the field interacts only with a group of slow ions \v z | <§C vt- When we turn the field off the resonance, the light interacts 
with faster excited ions. Collision frequency of those ions with ions in the ground state v is smaller, therefore the 
width decreases. It would be more proper to call this effect decrease of the broadening instead of the narrowing. 

To observe the narrowing experimentally, the probe-field spectroscopy can be used. One should measure the near- 
resonant absorption or gain of a weak probe wave as a function of its detuning on the same or adjacent transition. 
Another possibility is recording the spectrum of spontaneous emission in the presence of a strong resonant continuous 
field. Let us estimate two alternative effects independent of the detuning, namely, Stark broadening and change of 
velocity due to the interaction with electrons. The Stark effect is quadratic in the Aril spectrum. Its value Jl], p. 139] 
is about 100 MHz, which is small compared with both the homogeneous width T and the diffusion one kvT\Jvj\ 'j. 
The velocity change owing to scattering on electrons is as small as the ratio of the ion and electron thermal velocities 
VTi/vTe ~ 10~ 2 pl p. 155]. The ion-ion collision frequency decreases no more than by a factor of 2-3 at the detuning 
about kvx, therefore, neither alternative broadening mechanism can compensate the effect in question. The measured 
diffusion broadening in low-temperature argon plasma at zero detuning is by a factor of 3-4 (see [Q, p. 199]). At the 
detuning SI ~ kvx the narrowing factor expected is more than 10%, so this effect seems observable, see Fig. 1. 
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LIST OF CAPTIONS 

Fig.l (a) Distribution p^(v z ) of j -level population in ion velocity v z within the first order of the perturbation 
theory in the field intensity. Different curves correspond to distinct values of detuning O (from left to right Q/kvr = 
0; 0.5; 1; 1.5; 2; 2.5). Parameters are assumed to be T = 10~ 2 kv T , Tj — 10~ 3 kv T , v = 10~ 5 kv T . 

(b) Half- width on half a maximum of contour (v z ) as a function of detuning Q/kvx normalized to unit magnitude. 
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